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Abstract 

Star-shaped bodies are an important nonconvex generalization of convex bodies (e.g., linear 
programming with violations). Here we present an efficient algorithm for sampling a given 
star-shaped body. The complexity of the algorithm grows polynomially in the dimension and 
inverse polynomially in the fraction of the volume taken up by the kernel of the star-shaped 
body. The analysis is based on a new isoperimetric inequality. Our main technical contribution 
is a tool for proving such inequalities when the domain is not convex. As a consequence, we 
obtain a polynomial algorithm for computing the volume of such a set as well. In contrast, 
linear optimization over star-shaped sets is NP-hard. 

1 Introduction 

Convexity has been a cornerstone of fundamental polynomial-time algorithms for continuous as 
well as discrete problems [GLS88 . The basic problems of optimization, integration and sampling 
in W 1 can be solved efficiently (to arbitrary approximation) for convex bodies given only by oracles. 
More precisely, 

• Optimization. Given a convex function / : W 1 —* R, a convex body K specified by a mem- 
bership oracle and a point in K, and e > 0, find a point x* £ K s.t. f(x*) < mmx f(x). 
This can be done using either the Ellipsoid algorithm [YN76, GLS88 , Vaidya's algorithm 
|Vai96] or the random walk approach [BV044 ILV0 6a . For important special cases such as 
linear programming, there are several alternative approaches. 

• Integration. Given an integrable logconcave function / : M ra — ► R+ as an oracle, and e > 0, 
find A s.t. (1 — e) j f < A < (1 + e) j f . This is done using a Monte Carlo algorithm based 
on sampling logconcave densities |DFK91[ IAK91| ILV06bl |LV06a . 

• Sampling. Any logconcave density can be sampled efficiently |LV0 6a . The sampling algo- 
rithm is based on a suitable random walk. 

For the above problems and related applications, both the algorithms and their analyses rely heavily 
on the assumption of convexity or its natural extension, logconcavity. For example, for optimization, 
all the known algorithms use the fact that a local optimum is a global optimum. Similarly, a key 
step in the analysis of sampling algorithms is the derivation of isoperimetric inequalities, which 
are currently known for logconcave functions. Even the proofs of these inequalities (more on this 
presently) are based on techniques that fundamentally assume convexity. The main motivation of 
this paper is the following: for what nonconvex bodies/distributions, can the above basic problems 
be solved efficiently? 

In this paper, we consider a well-studied generalization of convex bodies called star-shaped 
bodies. Star-shaped sets come up naturally in many fields, including computational geometry 
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[PS85 , integral geometry, mixed integer programming, etc. |Cox73] . A star-shaped set has at 
least one point such that every line through the point has a convex intersection with the set. 
Alternatively, star-shaped sets can be viewed as the union of convex sets, with all the convex sets 
having a nonempty intersection. The subset of points that can "see" the full set is called the kernel 
of the star-shaped set. 

A compelling example of a star-shaped set is the "fe-out-of-m-inequalities" set, i.e., the set of 
points that satisfy at least k out of a given set of m linear inequalities, with the assumption that 
there is a feasible solution to all m. In this case the kernel is the intersection of all m inequalities. 
Another interesting special case is that of "fc-out-of-m-polytopes" , where we have m polytopes with 
a nonempty intersection and feasible points are required to lie in at least k of the m polytopes. 
These and other special cases have been studied and applied extensively in operations research 
[RP94, Mat94, Cha05]. Not surprisingly, linear optimization over even these special cases is -/VP- 
hard |LSN07| . 

This might suggest that the problems of sampling and integration are also intractable over star- 
shaped bodies. Indeed convex optimization is reducible to sampling. Our main result (Theorem [3]) 
is that, to the contrary, star-shaped bodies can be sampled efficiently, with the complexity growing 
as a polynomial in n, l/e,lnZ) and I/77, where n is the dimension, e is an error parameter denoting 
distance to the true uniform distribution, D is the diameter of the body and 77 is the fraction of 
the volume taken up by the kernel; we assume that we are given membership oracles for S as 
well as for its kernel K and a point xo so that the unit ball around xq is contained in S. (For 
the particular cases considered above, these oracles are readily available). The sampling algorithm 
leads to an efficient algorithm for computing the volume of such a set as well. We note here that 
linear optimization remains NP-hard even when the kernel takes up most of the volume. 

A reader familiar with sampling algorithms for convex bodies will recall that such an analysis 
crucially uses isoperimetric inequalities. Here we prove isoperimetric inequalities for star-shaped 
sets (Theorems [T] [2b. The key technical contribution of this paper is the proof of these inequalities 
and a new tool we develop for this purpose, which is also of independent interest. We refer to 
this tool as a thin decomposition of a set. The other crucial ingredients for efficient sampling 
(local conductance, coupling, etc..) extend naturally from the convex case to the star-shaped one. 
Therefore building on this new isoperimetry, we are able to show that the ball walk provides an 
efficient sampler for star-shaped bodies. 

In the rest of this section, we give some context for thin partitions. 

The common ingredient of most proofs of isoperimetric inequalities for convex bodies is the 
localization lemma, introduced by Lovasz and Simonovits |LS93] . The approach is based on proof 
by contradiction. If a certain target inequality is false in IR n , then there exists an essentially 
one-dimensional object over which it is still false. The proof is then completed by proving a one- 
dimensional inequality. This approach has been quite successful for convex bodies and logconcave 
functions and for proving many other inequalities in convex geometry. These, in turn, have played 
an essential role in the analysis of algorithms for convex bodies. 

However, this approach does not seem to work for nonconvex sets, since the resulting one- 
dimensional versions could be nonconvex or nonlogconcave (e.g., for star-shaped bodies, convexity 
holds along lines that intersect the kernel but is not required along lines that do not intersect the 
kernel). To overcome this, we use partitions of M. n induced by hyperplanes where each part is "long" 
in at most one direction. The overall proof strategy in applying the partition is proof by induction: 
we combine inequalities on all the parts to derive an inequality for the full set. The advantage of 
this (as opposed to proof by contradiction) is that a suitably strong inequality does not need to 
hold for every part; it suffices to hold for most parts. 
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1.1 Preliminaries 



Let S C 1" be a compact body. Define the kernel of S as K$ := {x : x £ S s.t. Vy G S [x, y] C S}. 
We say S 1 is star-shaped if Kg is nonempty and let rj(S) = vo\(Ks)/vo\(S). 

We denote the n-dimensional ball of radius r centered around a point x as B n (x, r). The ball walk 
with step size 5 in a set S is the following Markov process: At a point x in S", we pick a uniform 
random point in B n (x, 5) and move to the chosen point if it is in S and otherwise stay put. Let its 
denote the uniform measure on S and let a m denote the measure after m ball walk steps. For two 
probability distributions a, r, the total variation distance is 

d tv (o-,T) = supa(ct(A) - t(A)) 



1.2 Results 

We begin with two isoperimetric inequalities for star-shaped bodies, one parametrized using the 
diameter and the other using the second moment. 

Theorem 1. Let S be a star-shaped body with diameter D and rj{S) > 0. Then for any measurable 
partition (S%, S3, S2) of S, we have that 

vol(5 3 ) > ^d(Si,5 2 )min{vol(5i),vol(5 2 )} 

where d(Si,S2) is the minimum distance between a point in S\ and a point in ,S 2 . 

The above theorem is nearly the best possible as shown by a construction in Theorem [9] 

Theorem 2. Let S be a star-shaped body with rj(S) > and Ms = Es[\\X — fJ>s\\ 2 ] where [is is the 
centroid of S. Then for any measurable partition (S\, S3, S2) of S, we have that either 

vol(S 3 ) > ^vol(S) 

or 

vol(5 3 ) > -^^d(5i,5 2 )mm{vol(5i),vol( < S 2 )} 
16 V M s 

where d(S'i,S' 2 ) is the minimum distance between a point in S\ and a point in »S 2 . 

Next, we turn to the complexity of sampling. We assume that we have an oracle for the star- 
shaped body S, a lower bound on r] and an M-warm start o~o for the random walk, i.e. an initial 
distribution on S such that \/A C S, <jq{A) < Mtts(A). 

Theorem 3. Let S in M. n be a star-shaped body with kernel Ks and rj(S) > 77, and diameter D. 
Let its be uniform distribution over S and e > 0. Given a random point xq from a distribution o~q 
such that o"o is an M-warm start for tts, then there exists an absolute constant C such that, after 

Cn 2 D 2 M 2 , 2M 
m > ^ log 

steps of the ball walk with 8A ^^ -steps, we have dxvio'm^s) < e - 

Theorem 4. Let S in W 1 be a star-shaped body with kernel Ks and r]{S) > r). Suppose we are given 
membership oracles for Ks and S and a point xq with B n (xo, 1) ^ Ks Q S C B n (0, D). Then, for 
any e > 0, a nearly random point x from S can be produced using amortized 0*(n 3 /?7 4 e 2 ) oracle 
calls with the guarantee that the distribution a of x satisfies drv{, a ^s) < e - 

We note that up to the polynomial in r/, this matches the best-known bounds for sampling 
convex bodies. Due to page restrictions, many of the proofs appear in an appendix. 
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2 Thin Decompositions via Bisection 



Definition 5. Let S C M. n . We define S to be a compact body if S is compact, has non-empty 
interior, and satisfies cl(S°) = S, where cl(S°) denotes the closure of the interior of S. 

Let S C R n be a compact body. A decomposition of S is a finite collection V = {Pi, . . . , Pk} of 
compact bodies such that 

1. s = ufUPi 

2. Pi n Pj = dPi n dPj, 1 < i < j < k 

Furthermore, we define a decomposition V to be e-thin if each P £ V is contained in a cylinder of 
radius at most e. 

For completeness, we state without proof the following simple lemma. 
Lemma 6. Let S C W 1 be a compact body. 

1. Let N be a decomposition of S, and let S' C S be a compact body. Then 
N 1 = {cl((P n S)°) : P € N, P n S° + 0} is a valid decomposition of S' . 

2. Let N be a decomposition of S, and let N' be a decomposition of an element P € N . Then 
N U N' \ P is a valid decomposition of S. 

The following simple lemma from [LS93J will be used repeatedly. 

Lemma 7. Let f : W 1 — > M be integrable, n > 2. Then for any point z £ M n , and any 2- 
dimensional linear subspace S ofM. n , there exists a hyperplane H = {x : a T x = a T z}, with a G S, 
inducing halfspaces H + ,H~, such that it equipartitions f, i.e., 

/ f{x) dx= f{x) dx. 
JH+ JH- 

Theorem 8. For any integrable function f : W 1 — > K with supp(/) C S, S a compact body, and 
j fdx = ; and any e > 0, there exists an e-thin decomposition V of S such that each part P £ V 
is obtained by successive half space cuts from S and satisfies J p fdx = 0. 

Proof. Pick D such that S C M n (0,D). Since S is compact we know that D < oo. We start 
with the initial decomposition Vq = {S} of S. From this decomposition, we will inductively build 
decompositions Pi, . . . ,P n -i with the following properties. For each i, < i < n — 1, we have that 
for all P e Pf 

1. P is obtained from S via successive half space cuts. 

2. J p fdx = 

3. 3S C ]R n , an i-dimensional linear subspace of M n such that the orthogonal projection of P 
into S is contained inside of cuboid of side length at most 8 = ^f= . 

Assuming the above properties, one can easily see that each part in P n -\ is contained inside a 
cylinder of radius v^f = e i an d hence P n -\ is an e-thin decomposition of S compatible with / as 
needed. Hence, we need only show how to perform the induction step. 

Take P G P{, Q < i < n — 2. By assumption, there exists an i-dimensional linear subspace T 
such that tvt(P), the orthogonal projection of P into T, is contained inside a cuboid of size length 
at most 5. Since i < n — 2, we may pick a 2 dimensional subspace S orthogonal to T. 
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Let Q = conv.hull(-P) and let lis denote the orthogonal projection map from W 1 onto S. Since 
P Q S C M n (0,D) and JB„(0,L>) is convex, we know that Q C M n (0,D). Therefore U S {Q) C 
B n (0, £>) nS => vol 2 (<3s) < 7rZ) 2 . Let JV = {Q}. We perform the following iterative cutting 
procedure on N. Take an element E £ N. If vol 2 (ns(i?)) < <5 2 /2 stop. Otherwise, letting fi 
denote the centroid of Hs(E), we have by lemma [7] that there exists H = \x : a t x = a*/i}, where 
a £ S, such that j EnH - fdx = j EnH+ fdx = 0. Let Ei = E f) H~, E 2 = E f] H + . Now set 
N <— N U {i£x,-E 2 } \ £7. Since we are cutting through the centroid of Hs(E), and Ils(E) is 
convex, by Grunbaum's theorem we know that vol 2 (n s CEi)), vol 2 (n s (£ 2 )) < (1 - i)vol 2 (n s (^)). 
Therefore, after a number of iterations depending only on D, we will have that every element E G N 
has vol 2 (n 5 (£)) < f . 

Claim 1. Let E G iV. There exists v G 5, ||f|| = 1, suc/i t/iat widths (£^) = sup^g^f^x — 
inf^gs w T x < 5. 

Proof. Assume not, then note that the diameter of Ils(I?) is at least 5. Let [u, v] be a diameter 
inducing chord in ILg(E). Let to, z be points on opposite sides of [u, v] such that their distances from 
the line l(u, v) are maximum. Then the sum of the distances from w, z to l(u, v) is at least 5 and 
therefore the area of the quadrilateral induced by these four points is at least 5 2 /2, a contradiction. 
Hence there exists a direction v such that widths (E) < 5. □ 

Note then that the orthogonal projection of E into the subspace spanned by v and T is contained 
inside a cuboid of size length at most 5 as needed. 

Hence iV is now a decomposition of Q = conv.hull(P), such that each element of E G N 
has i + 1 orthogonal 5-thin directions. To transform N into a decomposition of P, we let N' = 
{cl((E n P)°) : E £ N, E n P° / 0}. By adding iV' to V i+ i, we complete the induction step as 
needed to prove the theorem. □ 

3 Application to Nonconvex Isoperimetry 

The benefit of Theorem [8] is that it will allow us to derive isoperimetric inequalities for high- 
dimensional sets without requiring convexity along every line. We show an application to star- 
shaped bodies. To gain some intuition, it is useful to understand what the obstructions to isoperime- 
try in the star-shaped setting are, as well as to understand why star-shaped bodies have good 
isoperimetry at all. The following Theorem illustrates what a "canonical" bottleneck looks like in 
the star-shaped setting. 



Figure 1: Star-Shaped Gluing of 2 Truncated Cones 

Theorem 9. Let H n = {x : x G = 0} and let (H~,H£) denote the halfspaces induced by 

H n . There exists an absolute constant C > 0, such that for all r] > 0, there exists a sequence 
of symmetric star-shaped bodies S n G M n centered at such that for all n > N„, we have that 
v(Sn) = ^(v) an d 

( ^Inf 1 ) \ 

vol n _i(fl„ nS n )<Ci _ ) I min {vol n (£T n S n ), vol n (if + n S n )} 
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Proof. Our strategy here will be to reduce the above statement to one about 1 dimensional dis- 
tributions on the real line. For each 77, we will construct a candidate sequence S n of star-shaped 
bodies which are rotationally symmetric about the x\ axis. Then by analyzing the cross-sectional 
distributions of S n and Kg n along the xi-axis, we will explicitly construct a 1 dimensional asymp- 
totic densities f v , fjf to which the cross-sectional distributions of S n and K$ n respectively converge. 
We will then establish the required isoperimetry and kernel volume constraints for the sequence S n 
and K$ n by direct computation on f v , fff . 

The geometry of our constructions is simple. As shown in Figure 1 previously, we will take 
two n-dimensional rotational cones with variance 1 along the x\ axis, truncate them at their ends 
removing exactly an r\ fraction of their volume, and glue them together at the truncation sites. 
Choose l n such that (1 



2l n < y/n(n + 2). Now let 



y/n(n+2) 



) n = 77. Since l n — > hi(^) we may choose N v such that for n > N v 



S, 



x : 



From here one can easily verify that the kernel of S n is 

In. X 



xf < 1 , 

^ =2 l ~ v^T2 



Xi G [0, l n ] 



K Sn 



L + x 



'■■ \ g*< £ 1 - ^m' x ' € 1 -'- 01 } U |- ^ g*?s 1 - t=i> 



xi G [0, Z n ] 



Next, a simple computation reveals that the cross-sectional distribution of S* n is 

x g [-l n ,0] 
x g [0, l n ] 



fn(x) 



' j y/n{n+2) , _ ln+x , n-1 _ 
2(1-7/) " V ^/n(n+2) ^ 

^/n(n+2) , _ ; n - x x n _i . 



1 



2(1-77) 

: 



otherwise 



Another computation, shows us that the cross-sectional density of K$ n relative to S n (we normalize 
by the volume of S n ) is 



1 y/n{n+2) , _ __^-x\n-l 
2(1-7?) n 11 ^^+2) J 
1 \/n(n+2) / 



2(1-7?)" 
: 



'(1 



Jn+x_\n-l 



V / «(™+2) 



X G [-Z„,0] 

X G [0, In] 

otherwise 



From here, one can easily verify that the sequence f n converges pointwise to the density function 
f„ : R -> R+ where 

'2(1^ : ^[-ln(I),0] 

AO*) = j w^rf x ■ x€[o,h^)] 

: otherwise 
Similarly the sequence converges pointwise to fff where 

xG [-ln(|),0] 
x€ [0,ln(i)] 



2(1-7?) 
77 

2(l-r?) 

: 



otherwise 
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Notice that < f v and that fff is log-concave. We get that 



/ f?(x)dx 



ln ^) 77 

~ X dx = 2 n . ' (1-7/) =7/ 



2(l-r/)7o 2(1-7?) 



The above computation shows that the volume fraction of the asymptotic kernel is indeed rj as 
required. Clearly the length of the support of f v is 21n(^). So now we see that the isoperimetric 
coefficient for f v is bounded by 



mm 



{j ^ f v ( X )dx, / °° /„(x)dx} " I " 1 - 7? " (1 - 7 ? )|SU P P(/,)| 



where |supp(/^)| denote the length of the support of /„. Clearly the isoperimetry computed above 
corresponds asymptotically to that of the partitions H~ P\S n , H^f]S n . To conclude the argument we 
need only justify Diam(5 n ) — ► |supp(/ r? )|. As is this is not the case, but this can easily be achieved 
by scaling S n orthogonally to the x\ axis by a factor of \. By doing this, we are collapsing the 
sequence S n onto the x\ axis, without changing the cross sectional distribution along the axis, and 
hence asymptotically Diam(5 n ) will converge to the to |supp(/) r? | as needed. □ 

From the above theorem and illustration, we see that the isoperimetry of star-shaped bodies can 
be strictly worse than in the convex setting where the isoperimetric coefficient is always f2 ( Diam(g) ) • 
In particular, from Figure 1, we observe how contrary to the convex setting we can get a V-shaped 
break in logconcavity of the cross-sectional volume distribution of a star-shaped body. On the 
other hand, as we will see later via Lemma [2~Tj the severity of these breaks is strictly controlled 
by the kernel of S. For reference, in Lemma |21| we show that the cross-sectional distributions of a 
star-shaped body satisfy a form of restricted logconcavity with respect to the kernel. The rest of 
this section will be devoted to proving isoperimetric inequalities for star-shaped sets. In particular, 
in Theorem [T] we show isoperimetry for star-shaped bodies in terms of the diameter and 7/ which 

in light of Theorem 



is optimal within a factor of O f i-^ ■ 



The next lemma forms the technical core of the isoperimetry proofs for star-shaped sets. Infor- 
mally, we prove that for any thin enough hyperplane cut decomposition of a star-shaped set S, the 
parts of the decomposition that intersect the kernel of S are "almost" convex. This will in essence 
allow us to apply the isoperimetric inequalities developed for convex sets to the "almost" convex 
pieces from which we will extract the isoperimetric properties of S. 

Lemma 10. Let S be a star-shaped body with r)(S) > 0, and let (Si, S3, S2) denote a measurable 
partition of S where vol (Si), vol (S2) > 0. Then for every e > 0, there exists a decomposition V of 
S such that 

1. VP G V, mo\(Si n P)vol(5* 2 ) - vol(5i)vol(5 2 n P) = 0. 

2. 3N CIV such that 

(«) ZreN > (1 " *MS) 

(b) VP G N, P is e-convex, i.e., there exists P' C W 1 a convex body, such that vol(P A P') < 
emin{vol(P),vol(P')}. 

First we state and prove a few technical lemmas needed to prove this. 
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Lemma 11. Let K\,K.2 Q M n be compact bodies and let vti,7T2 denote the uniform measures on 
K\ , K2 respectively. Then 

vo\(Ki A K2) < emin {vol(i^i), vol(K2)} =^ dt v (^1,^2) < e 

Proof. By symmetry, we may assume that vo\(Ki) < vol(i"C 2 )- Let /i,/2 : K n — > K + denote the 
associated densities with respect to m, 1^2- The subsequent computation proves the result: 

dto(7ri,7r 2 ) = r / |/i0) - / 2 (x)|cZx 

- 2 I vd(i^ + ( 1 ^ vol(^)vol(K 2 ) } 

1 / vol(^i)e evoKiTx) 

-H^(^ +V0l( ^ l) ^ 1 ? 

= i(e + e) = e 



□ 

Lemma 12. For S 1 C M n , K$ is a convex set. Furthermore, if S is compact then K$ is compact. 

Proof. If Ks = we are done. Therefore assume Ks 7^ 0, and pick x, y G Ks- Now take z G [x, y\. 
We need to show that Vp G 5, C 5. Assume not, then there exists p £ S, q ^ S, such that 
g G Since x,y £ S we have that C S. Furthermore we see that q is in the interior 

of the triangle defined by x, y,p. Let l(x, q) denote the line through x, q. Since q is in the interior 
of conv.hull {x, y,p} we must have that l(x,q) intersects the segment [y,p\ in some point r. But 
now note that r G S, since r G and by construction [x,r] % S, a contradiction. This proves 

the statement. 

For the furthermore, we assume that S is compact. To show that Ks is compact, we need only 
show that Ks is closed. If x is a limit point of Ks, we have a sequence {xj}^ C .Kg converging 
to x. Now to take any point p G S. We see that [a?i,p] Q S for all £ > 1, and we note that the 
sequence of line segments [xi,p] converge to [x,p] as i — > 00. By compactness of S 1 , we have that 
the limit segment [x,p] is indeed contained in S. Since this holds for all p, we see that x G Ks are 
needed. □ 



Proof of Lemma 10 (Near Convex Decomposition) . First we will show that we can find subset K T S C 
Ks that takes up most and the kernel and that lies deep inside it, i.e. that K r s + M n (0,r) C Ks- 

Formally, let K r s = {x : B n (x,r) C Ks} where r > 0. Let K s denote the interior of Ks- We note 

1 

that K s = U^-KTg. By the continuity of measure, there exists a positive integer j, such that for 
eo = j, vol (if £°) > (1 — e)vo\(K s ). Since Ks is convex, we know that vo\(K s ) = vo\(Ks) and hence 

vol(S) " vol(S) " ( 

Let / : M n — > E be /(x) = vol(<S , 2)lsi(^) — ^o\(Si)ls 2 (x) where ls^,!^ are the indicator 
functions for Si and S2 respectively. We note that f s f = vol(S , 2)vol(S'i) — vol(5 , i)vol(5 2 ) = 0. By 
Theorem |8j for every ei > 0, there exists an ei-thin decomposition V ei of S such that each part 
P G V is obtained by successive half space cuts from S and f p fdx = 0. We note that the condition 
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Figure 2: P is e-convex 



jpfdx = immediately implies condition (1) for V tl . For the time being we will assume that 
e i < \ e o and determine its exact value later. 

Let TV = {P : P G V ei , P n K e s ° / 0}. Since V ei is a decomposition of S, we note that U PeN P ~2 



and hence 



J2 PeN vol(P) _ vol(U PeiV P) ^ vol(gg>) _ 



vol(S) vol(S) ~ vol(S) 

We will now show that for an appropriately chosen e\ every P G iV is e-convex. Our strategy 
is as follows. We analyze a minimal cylinder C of radius €\ containing P, which exists by our 
assumption on V ei . We will use the fact that P contains a point deep inside the kernel to show 
that a subcylinder C of C is fully contained inside S. Lastly we will show that P' = C n P is a 
convex body whose volume is at least a (1 — e) fraction of the volume of P. 

Take P G iV. Let C be the cylinder of radius at most ei and let L C C denote the axis of C. 
Without loss of generality, we may assume that L is a subset of the x\ axis, i.e. 

C=|x:xGM",a<x 1 <6,^xf <e?|. 

By assumption, we have that P Pi -fT^ / 0, so pick c£ Pfl if^ . Since ei < ^eo, there exists 
d G L such that ||c - d|| < ei < \e Q . Hence B n (c, e ) C K s => B n (d, f 1 ) C Kg. Let 5 = jeo- 
Without loss of generality, we may assume that d = 0. Furthermore, by choosing a, 6 minimal 
subject to containing P, we may assume there exist points v,w G P such that -ui = a and u> i = b. 
By possibly rotating C, we may assume that v = (a, r, 0, . . . , 0) where < r < ei. By assumption 
on d, we know that t = (0, —5, 0, . . . , 0) 6 Ks- Therefore the line segment [v, t] C S. By a simple 
computation, we see that [v, t] intersects the x\ axis at v' = (^ga, 0, . . . , 0). Since G Kg, we also 

have that [v' , 0] G S v* = (^q^a, 0, . . . ,0) € S. By symmetric reasoning with respect to w, we 

have that w* = (^b, 0, . . . , 0) G S. 
Now, consider the subcylinder 

Claim 2. C" C 5. 

Proof. Take x G C". By symmetry we may assume that x = (e, /, 0, . . . , 0) where < e < jEpf^& an d 
< / < €\. Now examine the line l(x,w*). A simple computation reveals that l(x,w*) intersects 
the X2 axis at the point x* = (0, Sb _^_ ei ^ e f, 0, . . .). Now we note that 

5b , 5b 5b 5b 

:/ < 77 77 , -T- Cl < ~ 77" 777777 e l = 7T e l = 6 



5b -{5 + eiK ~ 5b-{5 + ei)e ~ 6b - (6 + ei)|=f^6 ei& 
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Therefore by assumption on d we know that x* G K$- Since x G [x*,w*], we have that x G 5 as 
needed. □ 

Now define P' := P n C 
Claim 3. P' is convex. 

Proof. To see this note that P is obtained from 5 via halfspace cuts, i.e. P = Sf)™ =1 Hi where 
each Hi denotes a halfspace. Now we see that 



TO 

p' = p n c' = s n C P| Hi = C P| Hi 



1=1 1=1 

since C" C 5. Since the intersection of convex sets is convex, we have that P' is convex as 
needed. □ 

Now note that P A P' = P \ P' = P \ C' . We will now show that for an appropriate choice of 
e±, depending only on 5 and e, we have that PAP' < evol(P') which will prove that P is indeed 
e-convex. In fact, letting P + = P n {x : x\ > 0} , P' + = P PI {x : x\ > 0}, we will prove that 

vol(P + AP|) < evol(P|) 

By symmetry, the same inequality will follow for the x\ < side, and by summing up the two 
inequalities the result follows. 

Let S(t) = {x : x G P+,xi = t}, and s(t) = vol„_i(5(t)). Now let b' = f=f^6 and let t* = 
&rgmax b , <t<b s(t). We have that 

vol„(P+ \ C) = t ' s(t)dt < f b s(t*)dt = (b-b')s(t*) = (^-] bs(t*) 
Jb> Jb> \o + e 1 J 

Now by construction the section 5(0) C Kg- I claim that 5(0) C Kp + . Take x G 5(0) and y G P + . 
Since x G w e have that [x, y] C 5. Now P + = 5P)^ 1 iJj. Clearly x,y G P =>■ x,y G PTj, for 
1 < i < m. Furthermore, since each Hi is convex, we have that [x, y] C PTj. Therefore [x,y] C P + 
as needed. Choose a G [0, 1] such that (1 — a)0 + at* = b' . Since 5(0) C Kp + , we see that 

(1 - a)S(0) + aS{f) C 5(6') 

Therefore by the Brunn-Minkowski inequality, we have that 

s(b') > ((1 - a)vol n _i(S(0))^ + ovol n _i(S(**))^) 
> a"- 1 vol„_ 1 (5(r)) > ( S j^-Y 'sit*) 



n-l 



Since P' + is convex we note that conv. hull {0, 5(6')} C P' + and hence 

1 1 / A \ n 

vol n (P|) > vol„ (conv.hull {0,5(6')}) = -b's{b') > - ( — — J 6s(t*) 
Now by choosing ei small enough such that 

+ eij n \d + ei 
we get that vol(P + AP|)< evol„(P^) as needed. 

□ 
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Using the above lemma, we now prove Theorem [T] 

Proof of Theorem [7] (Diameter isoperimetry). Let (Si, S3, S2) be a measurable partition of S. With- 
out loss of generality we may assume that vol(Si) < vol(S2) a = V °\^ S „ 1 \ < 1. Note that 

vol(S) > vol(Si) + vol(S 2 ) = ^^vol(Si) vol(Si) < — ^— vol(S) 

a a + 1 

Let V € be the decomposition of S with respect to (Si, S3, S2) as defined in Lemma 10 with parameter 
e. Let N denote the set of e-convex needles. Let N + = |p : P £ AT, vol (Si nP)> §^ v ol(P)j 
and N- = N \ AT+ . Since N = N+ U N- and by assumption on N 

vol(P) 



pen y ' 



we must have that either 



PEN- y ' PEN+ y ' 

Assume first that (a) is true. We will show that Si and Si take up a small fraction of most 
partition parts and that consequently S3 must take up a large fraction of S. Take P £ A^ _ , and let 
Sf = Si fl P, S2 = S2 n P, S3 3 = S3 n P. By assumption on V e we know that vol(Sf) = avo^Sf). 
Therefore we have that 

vol(Sf) +vol(S 2 D ) = (f^ti)volOSf) < -vol(P) 

by assumption on N~. Since vol(Sf) + vol(S<f) + vol (S3 5 ) = vol(P), we must have that vol(Sl') > 
^vol(P). Therefore, we have that 

V0KS3) v voKsf) v I^)>In- e ws) 

vol (5) " ^ vol(S) " ^ 2vol(S) "4 (1 e)7][b) 

v ' pen- y ' p^n- v ; 

Since d ( Sl ^ s ^ < i ; this proves the theorem for case (a). 

Now assume that (a) is not true. Then we must have that (6) is true to satisfy our assumption 
on N. Now take P £ N + . Our strategy here will be to derive isoperimetry for P using the fact 
that P is e-convex. By approximating the measure of P by that of its convex approximation, we 
will derive an isoperimetric inequality for P with an additive error depending on e. Since in this 
case Si and S2 take up a lower bounded fraction of P, we will be able to transform the additive 
error into multiplicative error by making e sufficiently small. The statement will follow as a result. 

So let Q be a convex body such that vol(P AQ) < emin{vol(P), vol(Q)}. We may assume that 
Q C conv.hull(P), since otherwise Qnconv.hull(P) is a convex body and strictly closer to P. Next, 
since Diam(P) = Diam(conv.hull(P)) we have that Diam(Q) < Diam(P) < Diam(S) = D. 

Let ttp,7Tq denote the uniform measures on P, Q respectively. Let Sf = Sf n Q, S% = Sf n Q 
and = Q \ (5? U S$). Since (Sf, Sf, Sf ) partition P, we note that Sf = (Sf DQ)U(Q\P). 
Then we have that d(S® , S%) > d(Sf , Sf) > e£(Si, S2). By lemma [IT we know that dt v {irQ,iTp) < e 
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and so we get that 

ttqOS?) = n Q (S[) > tt p (S[) - e > tt p (S[) - 3e 
ttq(S^) = ttq(S p ) > tt p (S 2 p ) - e > tt p (S 2 p ) - 3e 
vtq(^) < ttp(S^) + 6 = n P (S® n P) + e 

< ir Q (S® n P) + 2e = 7Tq(S 3 p ) + 2e < tt p (S p ) + 3e 

Since Q is convex, using the isoperimetric inequality proved in [LS93] we have that 

^^3) > dwqT mm \ 7r( 2 ( . 5 iW( ^2); > — -f) — ^HW^flw) 

Now bringing the above inequalities together, we get that 

MS 3 P ) + 36 > d( y 2) min {n P (S[ ) - 3e, ^p(5 2 p ) - 3e} => 
vtp(5 3 p ) > d{Sl D S2) min {^(S p ) - 3e, vr P (S p ) - 3e} - 3e => 
*p(# ) > Mtf ) " 3e) - 3e 



since 7r P (5f ) < 7r P (5 P ). Now choose e < min f§^} where e > 0. Since P G iV+ 

we have that vol(Sf) > |^!Vol(P) ttp(S[ ) > 5^. Hence 3e < f ^ < f 7r P (5f ). A simple 
computation now gives us that 

vr P (i 3 J > (1 - e ) 



2 L> a + 1 
Now vol(S p ) = 7r P (S' P )vol(P), so we see that 

V0l(5 3 ) > Yl V0 '(* S 3 P ) 



P&N+ 



a + 

PG-/V+ 

V a + 1 2 
>(l-e )(l-e) r ^-d(S 1 ,S 2 ) V o\(S 1 ) 

= (1 - e )(l - e)^d(S 1 ,5 2 )min{vol(5 1 ),vol( < S 2 )} 



Finally, letting eo — ► yields the result. □ 

We prove Theorem [2] following a similar proof strategy as Theorem [T] We need the following 
lemma about second moments. 
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Lemma 13. Let fx, . . . ,f m : IR n — > M + 6e densities with associated random variables X±, . . . , X m 
and centroids fix, . . . , /i m respectively. Let g = Y^ILxPifi be a mixture of the fas with associated 
random variable Y and centroid fj,. Then we have that 



l<i<j'<ra 



E[[|y-/,[| 2 ] = J>E[||X. 
i=l 

Proof. Since g is a mixture, we see that 

„ m . m 

/i= / xg(x)dx = y^Pi / xfi{x)dx = y^pmi 

Next, note that ||Y — /i|| 2 = (Y — /j, Y — /i). The following computation yields the result: 
E [||Y - /^ll 2 ] = / (x — fj,,x — /j,) g{x)dx 

(x, ») g{x)dx - 2 (fx, /i) + (/x, /i) = V" Pi / (x, a;} fi{x)dx - (fi, fjt) 

i=i jRn 

jil „ m /mm 

vv / ((x,x) - (/j,i,tH))fi(x)dx+y2p i (fH 1 fH) - (Ew^Ew 

7=t * n 7=t \£i J=t 

m mm 

5>e 



m 



E^ E i 

m 

E^ E i 

m 

E^ E i 

i=i 

m 

E^ E i 





- ml 




- ml 


\\Xi 


- ml 


\\Xi 


- ml 


\\Xi 


- ml 



] + E^ ^) ~ E^ 2 ^> ^) + E 2 P*Pi Mi> 

l<i<j'<m 



i=l 
m 



i=l 



2 ] + E^ 1 ~ ^ ~ E 2pi Pi (m,^j) 

l<i<j'<m 



1=1 
m m 



2 ] + E E Pi ^ (m,m) - E 2 W Mj) 

l<i<j<m 



t=l i=l 



2 ] + E P*Pi <A*i> A«i> - 2 PiPj (m, N) + PiPj (NiN) 



l<i<j<.m 



2 ]+ E PiPillm - mI 



i=i 



l<i<j'<m 



□ 



Proof of Theorem^ (Second moment isoperimetry) . Let (Si, S3, S2) be the measurable partition of 

S. We may assume vol(Si) < vol^) and so a = < 1- Let V e , N, N + , N~ be defined as in 

the proof of Theorem [T] Again as in Theorem [l] we have the cases (a) and (b). If case (a) occurs, 
then by the proof of Theorem [T] we have that 

vol(5 3 ) > i(l-e)r/(5)vol(5) 
as needed. So we may assume assume that we are in case (6), i.e that 



Pe7V+ 



vol(S) ~ 2 
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Now for each P G V £ , let ttp denote the uniform measure on P, fip denote the centroid of P, and 
let Mp = Ep[\\X — Hp\\ 2 ]- Now we note that tts, the uniform measure on S, is a mixture of the 
7rps, i.e. 

^ vol(P) 
^ vo(6) 

Therefore by Lemma [T3| we have that 

M ST V ° I(P) M 4 >T vol(P)vol(Q) ||j; n ||2> vol(P) 
PeP f v 1 {P,Q}c-p e v ' PeN+ v ' 

Let V = ^Pejv+ vol(P). By assumption 1/ > ^(1 — e)r/(5)vol(5), and hence 

Let N* = |p : P G N + ,M P < Mgj. Since Ep e 7v+ ^^r^Mp is an average of positive 

numbers by Markov's inequality we must have that 



£ vol(P)>^>i(l-e)77(fl)vol(5) 



PgJV* 

Now take P G JV*. By assumption on N* C JV, there exists Q a convex body such that 
vol(P A Q) < emin{vol(P),vol(<5)}. In particular, by the construction of Lemma 10 we may 
assume that Q C P. Let Q = P \Q, and let ttq,ttq denote the uniform measures on Q,Q 
respectively. We now see that 

vol(O) vol(Q) 

* p = ^UPj + ^(pf Q 

As done previously above from Lemma [13] we readily see that 

M P > V0] ^M Q =}► V °j^Mp > M Q (1 + 6)Mp > M Q 
vol(P) vol(Qj 

As in the proof of Theorem[l) let = S[ nQ,5 2 ° = n Q, and 5^ = Q \ (5? U S®). Since 
Q is a convex set, using the isoperimetric inequality proved in |KLS95| we get that 

iqQ\ ^ ^(^?> ^2) ■ f /oQ\ ^e<2\l ^ d(Si,S2) . f , c q. /o^hl 
^ } " ^ Q ( 1 ),7rQ( 2 ^ 2 x /(l + £ )Mp mm l^ ( ^ )'^^ } ) 

" { sil + eWs ) 2 d{Sl ' S2) min ^(S 2 Q )} 
Using the same analysis as in Theorem [T] the above inequality gives us that 

MS P ) > ( gl'^if] ) 2 52) min ^ p{s ^ " 3e ' " 3£ } " 3e 

Now choose 

1 



e < min 



(it) (svl+^S) d(Sl ' 52) (a^l)'(ll) Va + 1 



a; 
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for any eo > 0. By the same analysis as in Theorem [T] we get that 

Using the fact that J2peN* y0 K^) — |(1 — e )v(S)yo\(S) we get that 
vol(5 3 ) > vol ( 5 3 P ) 



pcn* 
> (1 - eo) 



> 



> 



8(l + e)M 5< 
(l- eo )(l- e )f / 77(5)1 



(1 + 6)2 

(1 - € )(1 - e)i / V(SV- 



PeN* 



a 



(l-6 )(l-6)l / V (S)I \ 

(l + e )| \IQ^/M~s) 



d(S 1 ,S 2 )—-vo\(S) 
a + 1 



d(Si,S 2 )\rol(Si) 
d( ( S , i,5 2 )min{vol(5i) J vol(5 2 )} 



Finally, letting eo — > yields the result. □ 

4 Conductance and mixing time 
4.1 Local Conductance 

Ball walk on star-shaped bodies could potentially get stuck in points with very small local conduc- 
tance. Here we prove that most of the points in a star-shaped body have good local conductance. 
First, we extend a lemma from [KLS97 from convex bodies to star-shaped bodies which leads to 
the proof of good local conductance. The proof is essentially identical to the case of convex bodies. 

Lemma 14. Let L be a measurable subset of the surface of a star-shaped set S in W 1 and let 

Sl ■= {(x,y) : x e S,y £ S, \\x - y\\ <r,xyHL^(p} 

Then the 2n— dimensional measure of Sl is at most 



r" +1 V0l n -l(L) 



7Tn-l 



71+1 

Proof. By the sub-additivity property of this measure, we have that if L = L\ + L 2 , then, 

vol(5 Zl U S L2 ) < vol(5 Ll ) + vo\(S L2 ) 

Therefore, it is enough to prove this in the case when L is infinitesimally small. Under this 
assumption, the measure ^{Sl) is maximized when the surface of S is a hyperplane in a larger 
neighborhood of L. Then the required measure is at most 



dydx 



x&S L :d(x,L)<r J y£S L :\\x-y\\<r 



dydxdt 



t=0 Jx£S L :d(x,L)=t Jy£S L :\\x-y\\<r 

vo\(Cap r (t))dxdt 

t=0 Jx£S L :d(x,L)=t 
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where Cap r (t) refers to the intersection of a ball of radius r and a halfspace at distance t from the 
center of the ball. Hence, 



l*(S L ) < [ [ vo\{Ca Pr {t))dxdt 

Jt=0 Jx£Sr.:d(x,L)=t 



t=0 J x&S L :d(x,L)=t 

ivol(JB n _ivV 2 -t 2 )dtda 



'xeL Jt=o 

rr 

vo' 



l n _i(L) tvol(B„_ivV 2 -* 2 )dt 
Jt=o 

vol n _i(L) f 7 r n _i(v'r 2 -t 2 ) n - 1 tdt 



n + 1 

□ 

Recall that the local conductance of a point x £ 5 is defined as Z(z) = v °l(jM x ' r ) ng ) 

K v 1 VOl(B„(0,r)) 

Corollary 15. Suppose a star- shaped body S contains a unit ball with the origin being inside the 
kernel. Then the average local conductance A with respect to ball walk steps of radius r is at least 



r\ n 

A > 1 - 



Proof. Let S and L be as in lemma 14 Then, if we choose x uniformly from S and u uniformly 
from B n (0, r), the probability that [x, x + u] n L ^ (j> is at most 

rvol n _i(L) 
2 frivol (S) 

Using the whole surface of S to be L we obtain that 

rvol n _i(<9S) 



A > 1 



2^vol(S) 



The corollary follows once we lower bound the volume of S in terms of its surface area. This volume 
can be written as the sum of the volume of cones whose apex is at the origin. Now, since the origin 
is present within the kernel, it can see every point on the surface of S. Hence, for each such cone 
C, vol(C) > — vol n _i (./4.c) ; where Aq denotes the base area of the cone. Summing this up, we get 
that vol(S) > ^vol n _i(d,S). □ 

The following lemma is the main result of this section. 

Lemma 16. Let S denote a star-shaped body containing the unit ball with the origin present inside 
the kernel and S r be defined as 

S r := {x e S : l(x) > ^} 

Then, 

1. vol(5 r ) > (1 - 2rVn)vol(5) 

2. vo\{K Sr ) > (l-2ry/n)\/o\(K s ) 

3. S r is star-shaped. 
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Proof. Using Corollary 15 we get that 



I f Ta IT) 

(1 - l{x))dx > (1 - -^-)vol(S) 



vo\(S) Js ~ 2 

V \ IT) 

E(l-l{x)) < -\~ 

^ ,'vol(B n (x,r) DS) 1\ r- 

Pr — v v , . . > - < 2rv/n 

1 vol(B n (0,r)) " V 

Pr(/(x)> 3 ) > (l-2rVn) 



> (l-2rVn) 



4, 
vol(5 r ) 
vol(5) 



Applying the same argument as above to the kernel of S r , we obtain the second inequality. The 
final conclusion is obtained as a consequence of 2. □ 

4.2 Coupling 

In this section, we prove that the one-step distributions of points close to each other and having 
good local conductance overlap by a good fraction. The proof follows the case of convex bodies 
closely since we are considering only points of good local conductance. 

Lemma 17. Let S be a star-shaped body and let u,v £ S such that \u — v\ < l(u),l(v) > I. 
Then 

d TV (Pu,Pv) <l + t-l 



Proof of Lemma 17 ( Coupling lemma). We prove the inequality in the case when both B n (u, r) and 
B n (u, r) are contained within 5. If not, then the considered case gives an upper bound and hence, 
we are done. 



d TV {P u ,P v ) = \{f \P u {x) - P v {x)\dx) 

Z Jx£B n (u,r)UB n (v,r) 

= hP u (u) - P v (u)\ + \\P u {v) - P v (v)\ + \ f \P u (x) - P v { X )\dx 

1 1 1 J xm n (u,r)nm n (v,r)\u\v 



+\ f {P u {x)-P v {x))dx + \ f 

z J xm n (u,r)rM„(v,r)\u z J xt 



(P v (x) - P u (x))dx 

l (v ) r)nB Tl (M,r)\i> 



= l(l-l(u)) + l(l-l(v)) + 

+1 [ (P u (x) - P v {x))dx + \ f (P v (x) - P u {x))dx 

1 AeB„(u,r)nl„(t),r)\ii z ii6B„(t),r)nB„(ti,r)\i) 

= 1-1+ 2vo | (B 1 n(0;r)) (2voKB,(0, r)) - 2vol(Ce^))) 
< 1-l+t 

□ 

4.3 Conductance 

Now, we bound the s— conductance of the ball walk on a star-shaped body. 
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Lemma 18. Let S C W 1 be a star-shaped body with diameter D such that rj(S) = r] fraction of its 
volume is present in its kernel. Then there exists a ball walk radius r such that the s-conductance 
<J> S of ball-walk of radius r is at least 2 iSnD • 



Proof. Let the radius r of the ball walk step be s/4i/n. By Lemma 16, this gives us that 

vol(5 r ) > (1 - ^)vol(S) 
Further, the fraction of the volume of the kernel of S r is 

Now, let A U A be any partition of S into measurable sets with vol (A), vol (A) > s(vo\(S)). Define 
sets 

A 1 :={x£AnS r :P x (A) < ^} 
A 2 :={xeAnS r :P x (A) < 

lb 

A 3 := S r \AM2 

Now, suppose that vol(^4i) < VO g g ^ . Then the conductance 4> S (A, A) is at least 

1 dx = — — -— -^-j- ^-vo\(An S r \Ai) 



mm{vo\(A),vo\(A)} J x£AnSr \ Al 16 min{vol(,4), vol(A)} 16 



^ ■ J [A\ ,,T U ^(vol(^i)-vol(S\S r )) 
mm{vol(j4), vol(A)} 16 

mmjvol^), vol(A)} 16 3 2 
^ ■ j ul\ irlu l7(T vol ^-^ vo1 ^ 

min{vol(il), vo\{A) \ 16 3 2 

> ^—w\(S) 

~ min{vol(A),vol(A)}32 1 ' 



> 



1 

32 



and hence we are done. Therefore, we may assume that vo\(Ai) > and vo\(A 2 ) > 

Consider u £ A\ and v £ A 2 . Then, 

d T v(Pu, Pv)>l- P U {A) - P V {A) > 1 - - 

8 

Using Lemma |l7| (t=l/8) , we get \u — v\ > and hence, d(Ai, A 2 ) > Now, using Theorem 
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[T]on the partition A\, A%, A% of S r , we get that 



$s > ? 7^ / Px(A)d,X 

~ min{vol(A),vol(A)} J A xK ' 
> 11 vo\(A 3 ) 



> 
> 

> 
> 



2 16min{vol(A),vol(^)} 
1 r t (S r )d(A 1 ,A 2 ) min{vol(^i),vol(A 2 )} 
2 5 4D min{vol(A),vol(I)} 
1 ?/(5r(l - s/2)) min{vol(A),vol(A)} 
29 8^/n~D min{vol(A),vol(A)} 
1 s(l - s/2)?? 
212 
1 S77 
2PnD 



□ 



Using Theorem |2j one can derive the following bound by proceeding similarly as in the proof of 
the above lemma. 

Lemma 19. Let S C W 1 be a star-shaped body with diameter D such that rj(S) = r] fraction of its 
volume is present in its kernel. Then for ball walk radius r = s/Ay/n, for any partition A, A of A 
satisfying vo\(A), vol(^4) > s(vol(S)), the s-conductance of A satisfies 

, r a\ ^ V ■ f vol (5) S I 77 

4.4 Mixing time 



2 9 \min{vol(A),vol(A)}' 2 7 n ]j M s 



Let its denote the uniform distribution over the star-shaped body. Let o m denote the distribution 
after m-steps of the ball walk on the star-shaped body. To relate the s-conductance to the mixing 
time, we use the following lemma from [LS93 . 

Lemma 20. Let < s < 1/2 and H s = sup^/^^ |<r (^4) — 7rs(^4)|. Then for every measurable 
A C W 1 and every m > 0, 

\a m (A)-7rs(A)\<H s + ^(l-^- 

s \ 2 

Proof, (of Theorem [3]) Suppose o"o be a starting distribution such that there exists M > 0, \/A C S, 
o~o(A) < Mtts(A). Now, by definition H s < M ■ s. Hence, using Lemma 20 and Lemma 18 

(2 2 \ m 

Replacing s by e/2M, for m > 2 ™ n ^ M2 log we have d T y{o m , tt s ) < e. □ 



5 Sampling algorithm 

To obtain a polynomial-time sampling algorithm we make the additional assumption that we are 
given an oracle to the kernel of the star-shaped body, a point xq in the kernel and parameters r, R 
such that B n (xo, r) lies in the kernel and the kernel is contained in a ball of radius R. The sampling 
algorithm proceeds as follows: 
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1. Use the algorithm of [LV07, LV06a] to find a transformation of the body S into isotropic 
position and obtain a random point xq in K$- 

2. Perform m ball- walk steps from xq on the transformed body S', for each desired random 
point. 



Clearly, by step 1 above, we have a ^-warm start for the ball-walk on S'. Now, by Lemma 
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to 

obtain a bound on the s-conductance, we need an upper bound on the mean square distance Mgi 
of the body S' . 

We next show that when the kernel is isotropic, the body is not far from isotropic. This will bound 
Ms 1 which along with Lemma [19] and Lemma 20 would prove Theorem [4] 

Lemma 21. Let S be a star-shaped body and let Kg be the kernel of S. For a vector v G M. n , 
\\v\\ = 1, define 

f s (t) = vol n _i({x : v T x = t,x G R n } n S) and f K {t) = vol n _i({x : v T x = t,x G R n } n K s ), 

the cross- sectional volumes for S and Ks in direction v. Then for x G supp(/^) and y G supp(/s) 
and a G [0, 1] we have that 

i i i 

fs(ax + (1 - a)?/)"- 1 > afK(x)™- 1 + (1 - ^fsiy)"- 1 

Furthermore, for all x, y G M. we get that 

f s {ax + (1 - a)y) > f K {x) a f s {y) l ~ a 

Proof. Let S(t),Kg(t) denote the cross-sections of S and K$ in direction v at t. Since x G 
supp(/is-),y G supp(/s) we have that Kg(x), S(y) / 0. Since Ks(x) is part of the kernel we 
have that 

aK s (x) + (1 - a)S(y) C S(ax + (1 - a)y) 

Therefore by the Brunn-Minkowski inequality we have that 

i i i 

afxix)™- 1 + (1 - a)/ s (y)«- 1 < vol n _i(aK s (x) + (1 - a)K s (y)) «-i 

i i 
< vo\ n -i{S(ax + (1 - a)y))*-i = fs{ax + (1 - a)y)^ 



For the furthermore, we note that the statement is trivial if either /k{ x ) = or fs{y) = 0. 
Therefore, we may assume that x G supp(/x),y G supp(/5). Since the harmonic average is always 
smaller than the geometric average, the statement follows directly from our first inequality. □ 

Lemma 22. Let S C W 1 be a star-shaped body with an isotropic kernel Ks such that rj = 
vo\(Ks) /vol(jS') . Then, in any direction v, for a random point X from Ks, we have 

E((^) 2 ) < 33 ^. 

Proof. Let v = (1,0, . . . , 0) T w.l.o.g. Consider the cross-sectional density fx induced by the kernel 
along v. Since Ks is isotropic, we have that /x(0) > g |LV06b] . 

Next let / be the cross-sectional density of the body S along v. It follows that 

/(o) > vMo) > 
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Let a = sup jx : f(x) < - , x < j and b = inf jx : /( 



b-a< 



vf K (Q) 

OLIJJ "> O/ . J \ U> I 

Suppose not, then 



■'•) ,mj y^,x > oj. We claim that 



t?/x(0) 



>?/x(o) 



/ 



V(x)dx>^%-a)>l. 



Now consider a point x = tb for f > 1. Then by Lemma 21 we have that 



f(b) = f{(l-l}o+ > ( ?? /x(0)) 1 -i/(x)^ =}► 

fiVfinMo)) 1 '* > /(*) ^(o) (^y)* ^ ^) 

The same inequality as above can be derived starting from any 6' > 6, and since for every such b' 
we have that /(&') < by continuity we have that for i > 1 



/(*) < ^(0) 



In 2t 



By a symmetric argument, the same bound holds for x = ta. Let p = J f(x)dx. The following 
calculation gives the result: 



E((i/x 



/ x 2 f(x) dx = p ( - [ x 2 f(x) J + A x 2 f(x)dx + / x 2 f(x)dx 

i ~oo \P Ja J J-oo Jb 



< 



pmaxja ,6 } + r nf K {$) 



a /■oo 

2„-ln2(x./a) dx+ / ^ In 2 (*/&) ^ 

/& 

+ 



x e 



pmax{a 2 ,6 2 } + r?/x(0)(a 3 + 6 3 ) 



+ 



2 In 2 (In 2) 2 (ln2) 3 



< (a 2 + b 2 ) + (2a 2 + 2o 2 )6 = 13(a 2 + b 2 ) < 13 



< 



3328 



V 2 fK(0) 2 ~ V 2 



□ 



Proceeding similarly as in the proof of Theorem |3j one can derive the proof of Theorem |4j 
Proof of Theorem^ (Polynomial time amortized sampling). Lemma 22 gives an upper bound on 



Ms < ^2^- Using Lemma 
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we get that the s-conductance of the ball walk on a star-shaped 



body S C M n with the kernel in isotropic position and rj = vol (Ks) /vol (S) satisfies 



r] 5 / 2 s 
2 21 n 3 / 2 



By the sampling algorithm of |LV07| lLV06a] . Step 1 of the sampling algorithm takes 0*(n 4 ) oracle 
queries. Since in step 2 of the algorithm, we started the ball-walk on S by choosing a random point 
from the kernel, and the kernel takes up at least an r\ fraction of the volume of S, a random point 
from it provides an (l/77)-warm start. Proceeding similarly as in the proof of Theorem [3] we get 
that after m > log ^ ball walk steps, dTv(o~m, ^s) < £• Hence, by performing m-steps of the 



ball-walk for each desired sample, we obtain the desired amortized bound. 



□ 
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6 Discussion 



We have presented isoperimetric inequalities and efficient sampling algorithms for star-shaped bod- 
ies, based on a new technique called thin partitions. Linear optimization is NP-hard on these 
bodies, even when the kernel takes up a constant fraction of the body. 

Theorem 23. Given a star-shaped polytope S, it is NP-hard to optimize a linear function over this 
body for any rj(S) < 1, even if S is well-rounded. 

This result follows easily from a theorem of Luedtke et al. |LSN07| and we include a proof in 
the appendix for completeness. Thus, quite unlike convex bodies, linear optimization is NP-hard 
over star-shaped bodies, but sampling remains tractable. Given the sampling algorithm, we can 
estimate the volume as follows: given an oracle for the kernel, we can sample from K$ and obtain 
the volume estimate for K$ using |LV07| lLV06a| : further, given that f](S) > rj we can also estimate 
i](S) using 0(^2) samples and output the product of the two as the estimate for volume. 

We believe the thin partition approach should be broadly applicable to proving inequalities in 
convex geometry, especially for inequalities that do not seem reducible to one-dimensional versions 
(e.g., the KLS hyperplane conjecture [KLS95 ). 
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7 Appendix 

7.1 Optimization over star-shaped body 

Here we prove that optimization over a star-shaped body is NP-hard. In particular, we reduce the 
clique problem to linear optimization over a star-shaped polyhedron. 

Definition 24. An instance of CLIQUE(k) is given by a graph G(V,E). The problem is to decide 
if there exists clique of size greater than k. 

It is well-known that CLIQUE(/c) is NP-hard. 

Note that the NP-hardness of optimization over a star-shaped body does not depend on the 
fraction of volume of the kernel. 



Proof of Theorem 23 We reduce solving CLIQUE(A;) to minimizing a linear function over a star- 
shaped body. Given a CLIQUE(/c) instance G(V,E), define variables x G 1". For each edge 
e = (i,j), define ip e G W 1 , 

1 if I = i or I = j, 
otherwise. 



Iff 



For every edge e, denote the set of constraints given by x > ip e as a block constraint. Consider the 
following formulation: 

-A fk\ 
Minimize/ (a;) = 2_j Xi ' satisfying a t least I I block constraints among: 

i=i ^ ' 

Me G E, x > ip e (1) 

Define the feasible polyhedron as S. 

Claim 4. The feasible polyhedron S defined by the above formulation is star-shaped. 

Proof. First note that any subset of block constraints among the given constraints define a convex 
body. Thus, the feasible polyhedron is a union of convex bodies. Further, x = (1 1 . . . 1) T satisfies 
all the constraints and hence, we have a non-empty kernel. □ 
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Claim 5. By adding new constraints, a new feasible star-shaped polyhedron S' can be created such 
that rj(S') is a constant. 

Proof. Clearly x% > 1 Vi G {!,..., n} is a feasible convex region contained in K$. Therefore, by 
adding constraints Xi < a V i G {1, n}, for appropriately chosen value of a(> 1), one can make 
rj(S') a constant. Note that the set 1 < ii < a V i £ {1, ...,n} is still a feasible convex region 
contained in K' s . Specifically, one can choose a = n, to see that 

>(—)">- 

\ n J e 

□ 

Claim 6. There exists a clique of size k in G(V,E), if and only if there exists x G S such that 
f(x)<k. 

Proof. Suppose the graph has a clique C{V',E') of size k. Then, consider x* G 1" such that 
x* v = 1 V v G V. Now, for every edge e = [i, j) G x* > i\) e is satisfied since, x* = ifif = 1 and 
Xj = ipj = l and x£ > for S V, ^ Since C is a clique, = ( 2 ) and therefore, ( 2 ) 
block constraints will be satisfied which implies that x* G S. It is straightforward to check that 

Suppose there exists x G S such that /(x) < k. The objective function fix) can be rewritten as 
™ n Fc_B-|F|>( fe )'{Si=i max e6F{V ; |}}- Hence, there exists FOE, \F\ > ( 2 ), such that the edges in 
F cover at most k vertices. Clearly, this is possible only when F defines a clique of size k. □ 

Suppose there exists an algorithm A to optimize over a star-shaped body P given as an oracle 
such that r](P) > c. Now, given an instance of CLIQUE(fc), we formulate the linear programming 
problem as above. Following claim [5] we can find an appropriate value of a and add constraints 
such that Tj(S) > c. Further, it is easy to make S contain a unit ball based on the value of a. 
Finally, the oracle queries can be answered by checking the number of block constraints satisfied 
by the point x. Hence, we may use A to minimize f(x). Let z be the objective value obtained 
by optimizing using A. Using claim [6j it is clear that if z < k, CLIQUE(/c) is a "Yes" instance, 
otherwise CLIQUE(A;) is a "No" instance. □ 
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